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in a Haniiltonian formulation. The blocks are two-dimensional. Choosing couplings 
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thus conserving solvability. The result can be cast again in covariant form. By way of 
example we apply this exact solution to calculate the evolution of entanglement. 
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1. Introduction 

The Dirac oscillator [1] has attracted considerable attention both due to its simple 
formulation and its analytical solutions. Yet we may ask, if we can go beyond the single 
particle problem and conserve some of that elegance and simplicity. A plausible scenario 
in particle physics is that of a Dirac particle bound by an oscillator and interacting with 
an external (free) non-abelian field of mesons modeled as particles with finite radius 
[21 13]. Specifically, if the internal group of the interaction is chosen to be SU{2) and the 
non-local potential dependence is linear in the position and momentum operators, the 
structure which leads to integrability is preserved. This linearity can also be viewed as 
an approximation to a more complicated field. The central subject of this paper will 
thus be to develop this scenario for one, two and three dimensional Dirac oscillators. 

To achieve this we shall reformulate Moshinsky's approach in a matrix 
representation, where we shall see that a simple transformation will take the Hamilonian 
formulation of the Dirac oscillator to a block- diagonal form with 2x2 blocks. This 
representation is not manifestly covariant due to its Hamiltonian form. The equivalence 
which follows from the spectral theorem is confirmed by the fact that we recover the 
well-known solutions. The coupling to the field, which carries isospin, will necessarily 
lead to 4 X 4 blocks, which can still be solved exactly. The solvable models are obtained 
by choosing the coupling in such a way, that the blocks remain uncoupled. We shall show 
that several interesting models can be derived in this way, and that the two-dimensional 
case has a special role as compared to the one and three-dimensional cases. A covariant 
formulation of the different cases is also given. Once we obtain the exact eigen functions 
and eigenvalues we can study dynamics. 

Several cases of interest are studied by choosing particular values for the parameters 
of the model. For instance, we can have a coupling to the field including explicit 
dependence on 7-matrices (Pauli coupling). The conservation of isospin + fermionic 
charge in the dynamics can be produced by restricting this example. The case of 
vanishing coupling with the external field except for its static part is included as the 
simplest extension of the Dirac oscillator. The absence of explicit fermionic rest mass in 
the hamiltonian is also interesting in the sense of mass generation due to the external 
field. 

By way of example and as an extension to the fidelity studies presented for the Dirac 
oscillator in [10] we take advantage of the additional degree of freedom and consider 
the evolution of entanglement between the oscillator and the field or, equivalently, the 
decoherence of the oscillator due to the field. 

In section two we will review the Dirac oscilator in one, two and three dimensions in 
a matrix representation and show, that it is solved by casting it in block-diagonal form, 
the blocks being two-dimensional. This will not only serve as a reminder of previous 
work and fix our notation. It will also set the stage for the generalization given in the 
following section. There we will implement the coupling to the isospin field and show 
that this can be done such that the block size is increased only to four. This leaves us 
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again with solvable models. In section four we give the covariant formulation of these 
models and in section five an in-depth discussion of the one-dimensional case. To get a 
feeling of the solvable dynamics we have found, we calculate in section six the evolution 
of entanglement between the field and the Dirac oscillator. We shall end by drawing 
some conclusions and giving an outlook. 



2. The algebraic structure of 1,2 and 3 dimensional Dirac oscillators 

The Dirac oscillator [T] was originally proposed as a way to introduce a linear potential 
which preserved the solvability of the resulting Dirac equation. We review some of 
the known results of this system by using a new notation. The treatment to be presented 
is related to the algebraic structure of the Dirac oscillator hamiltonians in 1, 2 and 3 
dimensions. The relation of this structure with the solutions of the associated eigenvalue 
problem is discussed. This should be helpful when we consider solvable extensions of 
this system. Let us start with the three dimensional case, originally considered by 
Moshinsky and Szczepaniak. The hamiltonian is 



if = a ■ (p + il3r) + m(3 (1) 

where we have adopted units such that the rest mass and the frequency satisfy muj = 1 
and h = c = 1. We use a representation of the Dirac matrices given by 



Defining the rising and lowering operators which act on big and small components of 
spinors in the form 



S+ = ( ° ^^j=a^^ l2, S_ = (S+)t, Es = /3 (3) 

and using creation and anhilation operators a = r + ip, a} = r — ip for the oscillator 
variables, the hamiltonian ([T]) can be written in 4 x 4 matrix form as 




I w 

or in algebraic form as 



if = S+S-a + S_S-at + mE3, (5) 

where the spin of the particle is given by S = I2 C?) cr. In ([5]) we have made explicit 
the separation of the spin and the generators S-t, S3 of SU{2) acting on big and small 
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components of spinors (or particle and anti-particle solutions under the appropriate 
rotation). We shall refer to S = (S+ + S_,z(S_|_ — S_),S3) as *— spin. It is worth 
to mention that the projection of the *— spin operator S is not only related to quasi- 
relativistic and ultrarelativistic components of states. When combined with oscillator 
operators, the *— spin projection corresponds to states with positive and negative 
energies. This can be seen from the Foldy-Wouthuysen transformation of this problem 

EE]. 

The solutions of the stationary Schroedinger equation associated to (I5j) were 
originally obtained in [1] by writing 

(6) 

and recognizing that the second row yields il)2 = {E + m)^^{(T ■ a);/'!. This in turn can 
be replaced in the first row of ([6]) to give 

(m~2E+ [{m -E)/{m + E)] (cr ■ a) (cr ■ a)'^) ^/-i = 0. (7) 
All this can be derived in a simple way by squaring the hamiltonian in ([6]), i.e. 

/ (^■at)(^.a) + m2 U \ p2 / V'l \ ..^ 

1^ (^■a)(^-at)+m2 ; 1^ V^J / ^ ' 

which explains the solvability of the problem, resulting in two decoupled equations where 
finite and infinite degeneracies appear. For a discussion of supersymmetry related to ([6]), 
see [8]. We shall prefer, however, the hamiltonian in form ([5]), since it is equally useful 
in deriving the solutions of ([6]) and it will prove helpful in further extensions of this 
system. The dependence of H on ladder operators as in ([5]) shows that / = a"!" • a + ^Ea 
is an invariant operator which commutes with the squared total angular momentum 
J^. Therefore, the eigenstates of / with total angular momentum j, projection mj and 
radial oscillator number n can be used to evaluate the 2x2 blocks of H and obtain the 
stationary solutions. We adopt the notation |n, (/, s)j, ruj) for oscillator states \nlm) 
coupled to spin s, where j is the total angular momentum with projection m^, n is the 
radial oscillator number, / is the orbital angular momentum and m its projection on the 
z axis. We denote *— spinors by |±). A pair of states with angular momentum j and 
such that /| ) = (2n + j — 1)| ) is given by 

= |n, (j - 1/2, l/2)j,m,)|-), |02) = |n - 1, (j + 1/2, l/2)j, m,)|+). (9) 

Another pair of states with the same angular momentum j but with /| ) = {2n+j)\ ) 
is 

103) = \n, (j + 1/2, l/2)j, m,)|-), \<P,) = |n - 1, (j - 1/2, l/2)j, m,)|+). (10) 
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The 2x2 blocks of H obtained from these states can be evaluated. We denote their 
elements by H{j, 2n + j- = {(t>i\H\(f)j) for i,j = 1,2 and H{j, 2n+j)ij = {(f)i\H\(j)j) 
for i,j = 3,4. From ([9]) and (ITOl) one obtains the two following matrices 



-m 



leading to the well known energies E"^ = iv? + 2{n + j) and E'^ = + 2n. Infinite and 
finite degeneracies come from these two blocks respectively. In connection with these 
degeneracies, it should be pointed out that a study of the spectrum of the Dirac oscillator 
[7j reveals a symmetry Lie algebra given by a direct sum of so(4) (a compact part 
with finite-dimensional unitary irreps) and so(3, 1) (a non-compact part with infinite- 
dimensional unitary irreps). 

The discussion on the algebraic structure above can be implemented directly in 1 
and 2 spatial dimensions as well. Let us introduce the notation if^*^^ for the hamiltonian 
in d dimensions. Take x = ri,p = pi, a = ai for the one- dimensional problem. For two 
dimensions, consider the chiral ladder operators 

Ar = ai+ ia2, Al = ai- ia2 = (Ar)* (12) 

with the properties [Ar, Al] = [Ar,{Al)*] = 0,[Ar,A^j^] = [Al,A{] = 4. The 
hamiltonians for one and two dimensional Dirac oscillators in original form are 

H^^^ = ai{p + i(3x) + m(3, (13) 

with tti = — (Ji, [3 = 0-3 and 

= J2 ai{pi + il3n)+m[5, (14) 

i=l,2 

with «! = —(72, a2 = — (Ti, /3 = (Ts. These hamiltonians can be cast in algebraic form as 
H^^^ = a+a + a-a'' + ma3 (15) 

= ^^Ar + (t_AJj + m(T3 (16) 

which exhibit some similarities with Both of them have a 2 x 2 structure: The spin is 
absent in one spatial dimension and a± corresponds to *— spin, while in two dimensions 
(T3 generates the U{1) spin. The solvabihty can be viewed again as a consequence of the 
invariants 

/^') = a^a + ^^3 (17) 
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in one dimension, and 

= + las, Js = ArA^j, - AlA{ + ^as (18) 

in two dimensions. We must note, however, that the two dimensional case exhibits 
some pecuharities which do not appear in the other cases. The conservation of angular 
momentum J3 comes from the particular combination of cr and A^ in H'^'^\ together 
with the absence oi Al^ A\. This absence is also responsible for the infinite degeneracy 
of all levels. On the other hand, the three dimensional example is manifestly invariant 
under rotations due to its dependence on S ■ a and S ■ Its infinite degeneracy comes 
from the interplay between the infinitely degenerate operator (cr ■ a)(<T ■ a) "I" and *— spin, 
dividing {H'^^^y in two blocks as in (|8]). These observations will be used in the following 
section. 

3. Solvable extensions in one, two and three dimensions 

Now that we have characterized the solvability of these systems as the reduction 
of H^'^^ to 2 X 2 blocks, we may readily extend it to 4 x 4 blocks by introducing 
appropriate potentials. The corresponding Lorentz invariance of such procedure will 
be discussed later. Consider a hermitean operator of the form $(r,p) as the potential 
to be introduced in the total hamiltonian. One has H'^'^^ = Hq'^^ + $, with Hq"^^ 
given by the dimensional Dirac oscillator. On physical grounds, this corresponds 
to a bound fermion perturbed by a momentum-dependent potential. The scenario in 
which such type of potentials can arise is discussed with more detail in further sections. 
For now, let us assume that $ may come from either a scalar theory interacting with 
the bound fermion through Yukawa couplings, or as the scalar component of a vector 
potential introduced by the minimal coupling prescription. These two possibilities are 
not equivalent. We may introduce also an internal group for this field, for example 
the SU{2) associated to isospin in the Yukawa theory or as the gauge group of a non- 
abelian field. In either case we consider the fundamental representation (doublets). 
In the three dimensional case, let us impose full rotational invariance of this field in 
order to preserve integrability, i.e. [$, Jj] = 0. In general, we shall see that a wide class 
of operators $(r, p) allows to retain not only integrability, but also solvability of the 
resulting model. To simplify the discussion, we linearize the field in the variables r, p 
(this approximation will be discussed in section 4). We simply have 



<|. = $o + $iS-a+<l>lS-a"^ (19) 

where ^q, $1 are constant operators. We denote by T the vector of Pauli matrices in 
isospin space and the corresponding ladder operators by T± = |(Ti ± 2X2). Choosing 
$1, $0 properly, we may write 
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$ = (r+S-a + r_S-at + 7T3j (20) 

with 7 a real parameter. We may also write the more general 

$ = (A + Sa^) (T+S-a + r_S-at + 7r3) (21) 

parametrizing the two possibilities for $ mentioned above through the real constants 
A, B. For y4 = we have a field which corrects the rest mass of the fermion, while 
i? = represents the SU{2) potential = ($, 0, 0, 0). See section 4. 

Lower dimensional examples should follow the same pattern to preserve 
integrability. Thus, we propose the extensions 

= (T+a + (T_a^ + mcTg + {A + a^B) {T+a + T^a^ + 7T3) (22) 
H^^^ = a+An + cr_4 + ma^ + {A + a^B) (t+Ar + T_AJj + 7r3) (23) 



H^^^ = S+S ■ a + S_S ■ at + mSa + {A + ^^B) (t+S ■ a + T_S ■ a^ + 7r3) . (24) 

With these extensions, it is evident that the new invariants for one, two and three 
dimensions are 

I^'^ = a^a + ^a; + ^T^ (25) 

/(2)=A«4 + ia3, J3 + ^Ts = AnA^R-ALAi + ^a, + ^Ts (26) 

/^'^ = at.a+^S3 + iT3, J = atxa+S. (27) 

Before discussing the solutions of the Schroedinger equations associated to f l22]|23|24p . 
our new hamiltonians deserve some comments. It is important to note that the 
solvability of the resulting Schroedinger equation for these hamiltonians resides 
completely on the existence of I^'^^ (which in this case means superintegrability). For 
instance, in the three dimensional case one may consider any potential of the form 
$ = F(T+S • a + T_S ■ at + 7T3) where F admits a power expansion. Evidently, 
= 0. A suitable group of states can be used to evaluate the 4x4 blocks of H. 
We describe this procedure in subsection 3.1, restricting ourselves to the linear case for 
simplicity. 

Despite the fact that the three extensions fl22| [23l [2^ exhibit a similar structure, 
the two dimensional case is, as before, special. In our strive to preserve integrability, we 
have introduced a field $ which depends only on Ar and leaving Al and out 
of the game. The main consequence of this construction resides in the permanence of 
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the infinite degeneracy of tlie system (energies do not depend on A^^A^), as well as a 
modification of the conserved angular momentum J3 by the addition of isospin (second 
equality in f l26l) ). The special nature of the original H^'^^ is responsible for restricting the 
type of linear extensions and the resulting spectra. In contrast, the three dimensional 
extension breaks the infinite degeneracy through T± and conserves the angular momenta 
J^, J3 without modification. We discuss this through the eigenstates and energies of H^^^ 
in the following section. 



3.1. Solutions of the three dimensional case 

Now we turn to the problem of finding the eigenstates of H^^\ For this, we simply 
consider states which satisfy /| ) = ),J^\ ) = j(j + 1)| ). These are four 
in number. Then we proceed to evaluate the 4x4 matrix H{N,j) = ( \H^^^\ ). 
Introducing subscripts S and T for *— spinors and isospinors, we write such states as 



K) = |n,(j + l/2,l/2)j,m,)|-)E|-)r (28) 
|0^) = |n, (j-l/2,l/2)j,m,-)|-)s|+)T 
10^) = \n-l, {3 - 1/2, l/2)j,m,)|+)s|-)T 
K) = In - 1, (j + 1/2, l/2)j, m,)|+)s|+)r 

where, as before, n is the oscillator radial number, j is the total angular momentum 
and vfij its projection in the z axis. These are eigenstates of J*^^^ with eigenvalue 
iV = 2n+j-l/2. The resulting 4 x 4 blocks of if with elements H{N,j)ki = {(l)^\H\(j)f) 
are 



-m 



-m 



\ 



{ A-B)^ {A-B)p{n + j) -p{n + 3) 
{A - 5)7 


/2n 



{A- B)p{n + 3) 



An +3 






m - (A + 5)7 

{A + B)y/^ 



m 4 



\ 

72^ 

^B)y/2^ 

(A + 5)7 J 



(29) 



and the secular equation \H{N) — E\ = can be solved explicitly using the formula 
for the roots of a quartic polynomial. However, we shall contempt ourselves with the 
explicit solutions of the 1 + 1 dimensional example in further sections, as the blocks 
of H^^'^ and H^^^ are quite similar when j = 1/2. It is also evident that the infinite 
degeneracy is now broken, since one cannot reduce H{N) to smaller blocks where only 
n appears. The exception to this ocurrs when A = B = 0, which obviously recovers the 
usual Dirac oscillator. 



4. Lorentz covariant formulation 



Here we discuss the specific form of Lorentz covariant external fields producing the 
solvable extension we have established. Before proceeding, we should mention that the 
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potential $ in one and two dimensions indicated in f l22|[25]) . is spin independent. This 
makes our task easier, since $ can be directly related to the external field as shall be 
indicated. The three dimensional case, however, should be treated separately due to 
the presence of S in the potential $ considered in i^^. Therefore, we shall consider 
three alternatives for generating $: a parity violating coupling to an external non-local 
vector field (mediated by 75), a Pauli coupling to a non-local field tensor and a Yukawa 
coupling to a non-local scalar field. All of them shall have an internal symmetry group 
SU{2) as indicated before. The alternatives involving a non-local gauge field can be 
brought together when linear approximations of fields are considered, as we shall see. 



4-1- 3+1 dimensional case 

Here we use the metric tensor 77 = diag (1, — 1, — 1, — 1). The original covariant 
formulation of the Dirac oscillator [TT] considered a four vector u such that, in some 
inertial frame, = (1,0,0,0). For the Dirac oscillator of many particles, such vector 
was related to the center of mass of the system. In the present case, the meaning of u 
is given through the introduction of the anomalous (Pauli) coupling to an external field 
tensor. The equation leading to the Dirac oscillator hamiltonian ([1]) is 

[iM + ^IvU^'rl) + m\i\) = 0, (30) 
where r^ = — {r^Uy)u^ and 



7i 



ia^ 







7o 



(31) 



Eq.( l30|) can be cast as 



[7^p'^ + m + S^,FnV^ = (32) 

with F^^ = u^r^ — u^r^ and S^^p = i [7^,71/]. Thus, the meaning of u is obtained by 
means of the divergence of the Faraday tensor as d^F^^ = —u", i.e.u is an external 
current. It must be mentioned that another possibility to write a covariant form of the 
interaction comes through the dual tensor. The anomalous coupling with an electric 
field in the form "^^S^yF^^ with F^'^ = e'^'^'^^ F^p, produces the same result. 

With the aid of the vector we can introduce more interactions in a covariant 
way. As announced before, a non-local, non-abelian field tensor J-"'"^ = J2i=iTiJ-' ^'^ can 
be introduced in the equation by means of the Pauli coupling. The theory underlying 
this type of field can be found in many references [13], |13] and was motivated by 
the possibility of introducing high order derivatives in field actions. From all these 
possibilities, we choose the path followed in [I5] by considering that the field acts on 
wave functions as a bilocal kernel, i. e. it can be represented as an operator depending 
on r^,p^ (but not S). The integro-differential equations obeyed by F^'' are also in [15] 



Dynamics of a Dirac oscillator coupled to an external field 



10 



and can be used to elucidate the nature of the external source generating such field. 
We shall proceed to find the external current producing such fields after we give 
the expressions for the interactions. Let us now consider a linearization of the field 
in the p, r variables: an approximation which corresponds to the first corrections due 
to internal structure of particles described by the non-abelian field. This means that 
both wavelength and particle structure are small compared to the scale of the system, 
which is dictated by the oscillator frequency. The expansion of the field in r^,pj^ can be 
truncated by keeping linear terms only. We propose 



= ^''"^'u^r^p (33) 
= e^-^PuxPxp (34) 
= 0, (35) 

for which the Dirac equation reads 

i^^p^ + m + Sp^F^" + BS^.T^'']^ = 0, (36) 
This equation reduces to the hamiltonian form given in (!24l) with A = 7 = 0, as can 



be verified by setting = (1, 0, 0, 0). To produce a term 7E3T3 in the hamiltonian, we 
simply recognize that such a potential corrects the rest mass of the Dirac particle. A 
non-abelian scalar field can be introduced for this purpose in analogy with the one and 
two dimensional cases explained below. Notably, one could also consider the dual tensor 
with components 7^^" = e^^Q,^f" and an interaction of the form 'j^S^^^'^'^ in order to 
produce fl2il) . 

The other case, i.e. B = 0, A = 1, comes from considering again a parity violating 
coupling when we use the minimal coupling prescription. Such a form is encountered in 
field theories where chirality plays an important role [16], e.g. the electroweak model. 
Again, we need non-locality in our fields. By using a non-abelian 4-potential (e.g. a 
current vector with isospin 1/2) of the form Af" = ELi ^i^i, where 

A^,=r^^, A^,=p^^, A^ = 0, (37) 
we can see that the equation 

[^,p^ + m + S,,F^'' + ^757;.^^]^ = (38) 

produces the desired potential in the hamiltonian for 5 = 7 = 0. The term with 
T3 does not appear in this case, but can be generated by the unitary transformation 
if)' = exp{iTsu,yr'')ilj. This leaves p± invariant, while the matrices Ti,T2 in the 
hamiltonian simply rotate. 

The nature of such field can be elucidated by inserting our in the corresponding 
non-local field equations and finding the external current producing it. For instance. 



Dynamics of a Dirac oscillator coupled to an external field 11 

using ^^'^ one has ^^'^ = u^{r'^Ti + p'^T2) — ^, v and its expression in terms of a 
4-potential (see (12) in [15j) reads 

^'^^ = i{[p^, B^]-^i^v) + [B,, B,]. (39) 
This yields B^ up to gauge transformations as 

= u^{^r^r\Ti + r^p''^T2) Bihnear in p, r. (40) 
The field equations give the current as 

/=^[p„^n+[i^.,^n (41) 

= —u^Ti + + trilinear terms in p,r . 

So far, no approximations have been made. If we focus on the independent and linear 
terms we find a) a constant current as in the Dirac oscillator case, b) a spatial part 
given by the momentum of the particle, as if the imposed external current would 'follow' 
the electron. 

In order to put the cases A ^ 0, B together, we may consider a 4-potential 
given by the combination AA^j^ + BB^, with A^, B^ given above. The linearity in p and r 
invoked in our argumentation (e.g. in the case where only the first non-local corrections 
are considered and field intensities are small compared to the oscillator frequency) 
implies that the interaction 'y^'y'^{AA^ + BB^) can be approximated by ^757^^74^ when 
inserted in the Dirac oscillator equation. The anomalous coupling 'j^S^^^^ with the 
combination of A^ and i?^ yields the one produced by B^ when we ignore quadratic 
terms in p,r and in the constants A,B. We would like to stress that despite the fact 
that the linear approximation is helpful to put both cases together, one may consider 
each case separately without approximations in the light of (I39|40p . 

The physical origin of our potentials, as we have seen, depends strongly on 
the validity of non-local theories and the possible existence of extended particles as 
discussed by Yukawa (composite particles described effectively by such fields is another 
possibility) . 

4-2. 2+1 dimensional case 

Here we construct a field producing the two dimensional hamiltonian (1231) . Greek indices 
now run from to 2. Since spin does not appear in the resulting potential $, we identify 
it with a non-local scalar field, according to Yukawa's formulation in [2], [3]. In the 
presence of isospin T, the potential has an expansion of the form 
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$ = . T = ^ T, I d^k^'Un^'^fjik, l)(^^^-'(MPp\ (42) 
i=i 

where fc*^*-' are 2+1 dimensional Lorentz vectors representing the momentum of a scalar 
particle of mass M, while the vectors represent the transversal radius. We shall 
consider only one Fourier component. According to Yukawa [3], the corresponding 4- 
vectors obey 

= fc(2)fc(2)M = (43) 

= /(2)/(2V = (44) 
= = 0. (45) 

where A is the radius of the scalar particle. Again, linearization in p will be necessary 
and it is related to the first corrections to a scalar particle (or a doublet of particles, 
since we use SU{2) as isospin) with small A. Now, we can choose a frame of reference 
where the set of vectors above is given by 



A;« = {^{kWy + M\ k^'\0) = (0, 0, A) (46) 

= (y'(A;{2))2 + M2, 0, k^'^) = (0, A, 0). (47) 



and define 



1^" A J^^ V ^ 

Finally, we give the field $ = ■ T as 

0^ = «,(iVm _ /(ly ^ 02 = 4'V^-/iV, 03 = 7- (49) 

The equation 

[-f^,{p'' + i-fuu''r'^) + m + = 0, (50) 

with 7o = o"3,7i = icri,72 = —i<J2, yields the two dimensional hamiltonian in ( 123|) for 
A = 0, B = 1. The frame of reference where this is possible, is dictated by the unitary 
time- like vector u chosen as = k^^'' — k^^^ = /c^^^ — K,^fl'\ The case A = 1,5 = is 
obtained by considering a non-abelian gauge field whose scalar component is written in 
the same spirit, i.e. 

[if.ip'' + iluu'^r'^ + Af,) + m]ip = (51) 
A^ = u^<^ (52) 
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4-3. 1+1 dimensional case 

We proceed along the same line by introducing a non-local scalar field. Greek indices 
now run from to 1 and Dirac matrices are 7o = o"3, 7i = —i<y2- We use the momentum 
and radius vectors fc^, in a frame of reference such that 

A;^ = (m,0), /;. = (0,A), (53) 
define = (y-)^^ and set the field components as 



4>i = K^rf", ^2 = If.p'', h = 7- (54) 

Using the same Yukawa coupling and choosing = we obtain the desired potential 
in the hamiltonian for A = 0. The case i? = is obtained as before, i.e. by using a 
gauge potential = 



5. Detailed analysis of 1 + 1 dimensional extensions 

Now we concentrate on the 1 + 1 dimensional case for simplicity. Furthermore, we shall 
find a similarity between the block ( l29l) for j = 1/2 and the corresponding block for 
the one dimensional problem. Let us write the Dirac equation in hamiltonian form 
by choosing = (1,0). It is useful to introduce a coupling constant a for the field 
components 0i,02 in fl54|) . The resulting Schroedinger equation in terms of ladder 
operators reads 



a+a + cr_a^ + mas + {Aa^ + B) {aTj^a + aT_a^ + 7 



.dip 



(55) 



for which the usual Yukawa coupling is obtained when A = 1 and 5 = 0, while for 
A = 0, i? = 1 we obtain the zero component of the non-local potential. Now we proceed 
to find the solutions of the stationary version of ( 155|) by noting that the operator 



/ = ata+-(a3 + T3)-l (56) 

satisfies [if, /] = 0, i.e. J is a conserved quantity and it gives the set of states for which 
H is reduced to a block diagonal form. Such states are formed by harmonic oscillator 
states \n) and states |rcr) for isospin projection r = ± and *— spin projection a = ±. 
The states are given by 

|0n = |n + 2)|--) |0-) = |n + l)|-+) 

103) = 1^ + 1)1 + -) m) = \n)\ + +) (57) 
with n > and satisfy 
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m) = nm, 



1,...,4. 



(58) 



The Hamiltonian H in ( l55|) has a matrix representation 

f Ho . . . \ 
Hi ... 

H2 



H 



V 



where the diagonal blocks H^ have matrix elements {Hn)ij 
form 



(59) 



)7|i/|<^"). They have the 



Hn 



f-m-{B-A)-f a{B-A)^nT2 
a{B-A)y^nT2 -m + {B - A)-^ 





V 



VTr+2 




v/^rr2 
VnTT 

m-(5 + A)7 a{B + A)^/nTT 
a{B + A)y/n + 1 m + {A + B)-f 



(60) 



except for the lowest states forming a singlet and a triplet which we now indicate. 

The state |0, ) forms a singlet of /, z.e. the only state with eigenvalue equation 

/|0, ) = -2|0, ). Its energy is given by H\0, ) = -(m + {B - A)'y)\0, ) 

trivially. The states |1, ),|0,H — ), |0, — h) form a triplet of / with eigenvalue — 1, 

i.e. ) = — 11, ) (similarly for the other two). They give a block of the form 



H_ 



( -m - (5 - y4)7 
a{B - A) 

\ 1 



a{B - A) 
-m + {B - A)'y 




(61) 



1 \ 



m- {B + A)-f I 

leading to the lowest energy levels as solutions of a cubic polynomial. These two cases 
represent corrections to the vacuum energies, as can be seen by choosing vanishing 
interaction parameters A, B. In the following, we shall concentrate on general states 
described by fl57|) . ^.e. states satisfying fl58|) for values > 0. The spectrum can be 
obtained from the 4 dimensional secular equation \Hn — E\^\ = using ( l60l) . with I4 
the 4x4 identity matrix. This, of course, lies on the fact that the formulae for the roots 
of a fourth degree polynomial are known [20] . Eigenvectors can be written explicitly as 
well. 

For the sake of simplicity we study particular cases and discuss their physical 
importance. 



5.1. The case A = 0, B = 1 



This case corresponds to gauge fields acting on the Dirac particle. Conversely, one could 
put B = 0, A = 1, with analogous solutions except for the presence of some signs in 
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the parameters a, 7. The situation is then quite general from the point of view of the 
solutions. The resulting eigenvalues flM|) exhibit a dependence on n such that the factors 
in front of it get corrected by the presence of 0,7. The oscillator vacuum energy n = 
is modified as well. Therefore, both frequency and rest mass acquire corrections due to 
the interaction, as expected. The blocks of H are given by 



Hr, 



( -m - 7 a^/n + 2 + 2 

Vn + l 

m — 7 a\/n + 1 



a\/n + 2 —m + 7 




V 



^/n + 1 a\/n + 1 m + 7 / 



(62) 



The following definitions are helpful 



p = - [^2n + 1) (1 + a^) + 272 + 2m2 



g = 16m^ [yr? 



1 + 7')) + fl+47') +2fl + 4m2 



^7' ) + 



16 (1 + (o? - 2) + 7^ - (1 + 27^) + a^) (n + 2) + 
16 (1 - + a^) (n + 2f 
r = 108( ((n + 2) (1 - o?) + - 7^) ((n + 1) (l - a^) + 



1/3 



3 + y/r'^ — Aq'-^ J 
The eigenvalues can be expressed as: 



+ 



1 fr + ^/r^'^^Aq 



(63) 



Ei{n) 



/p±s , (-1)' 
2 ~'~ 2 



1,2 



(64) 



3,4 



We denote the eigenvectors of without normalization by such that -ffnl^ 
Ei{n)Vi. Their components Vji for j = 1, ...,4 can be written explicitly as 

Vu = aV^rT2 [{n + 1)(1 - a^) + (m - Ei)^ - 7^) 
^2* = (m + 7) ((n + 2)(1 - a^) + + - 7^) - 

((n + 2)(1 + a^) -a^ + {m + 7)^) + (7 - m)^^ + E'f 
V3^ = a((2n + l)Ei-m-7) 
V^4i = Un + 2){a^ - 1) + (7 + - 



m 



(65) 



5.^. T/ie case A = 0, Ba = 1, m = 7 

We put this case separately, although it can be reached from the solutions in the last 
subsection. For this set of parameters, it is evident that a new symmetry arises in the 
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form of permutational invariance cr -H- T and conservation of the total spin magnitude 
(<T + T)^. In the context of particle physics, such a conservation law indicates that 
isospin + fermion charge is a good quantum number. The characteristic polynomial is 
factored in qubic and linear expressions as a result of the antisymmetric singlet and the 
symmetric triplet of the total spin cr + T. In fact, the cubic part is a depressed cubic. 
The spin zero state is infinitely degenerate with zero energy. We use the states 



\X1) = + 1)(| - +) - I + -)) = |n + 1)|0, 0), Ix^) ^ \n + 2)\ - -) ^ \n + 2)|1, -1), 

IXs) ^-^\^+ -+) + ! + -)) = 1^ + 1)|1'0)' 1X4) = \n)\ + +) = \n)\l, 1) (66) 

where \s,ms) are states of the total spin cr + T with s = 0, 1 and its projection 
rUg — —1, 0, 1. It is straightforward to compute {Hn)ij — (x^|-H"|Xj)- We obtain 



Hr, 









27 







^2(^ + 1) 








(67) 



^J2{n + l) yJ^Mn^^ 

y ^2(n + 2) -27 ^ 

This yields eigenvalues which are similar to those for the previous case, but with 



5.3. The case a — Q 

As a simple application of our relativistic model, one can eliminate the presence of the 
field, except for its static part. For a perturbed confined particle, this case only takes 
into account a static correction to its rest mass and zero point energy. As we shall see, 
the oscillator number is not affected by factors involving 7, leaving the Dirac oscillator 
frequency untouched. One could state that the simplest case of mass generation (even 
at the level of classical fields) is thus achieved. The matrix becomes 



V 



-m 



-{B- 


\/n + 2 




Ah 





-m + {B 




Ah 



Vn + l 



Vri + 2 


m- {B + Ah 








Vn + l 


m+{A + Bh J 



with eigenvalues 



(68) 



Ei{n) = -S7 - + 2 + (m - A-ff, £2(71) = -S7 + ^n + 2 + {m- Aj)^ 
E^in) = 57 - Jn + l + {m + A^y, E^in) = B-f + Jn + 1 + {m + A-f)^. (69) 
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The eigenvectors of Hn are given by the rows of the following matrix 



/ Eiin)+{A+B)'y-m 

E2in)+{A+B)-y-m 
Vn+2 I 





V 










E3{n)-{A+B)-f-m 
/n+1 



E4{n)-{A+B)'Y-m 
Vn+1 



1 
1 

7 
7 



\ 


1 

1 ; 



(70) 



5.4- The case m = 7 = 



This corresponds to a confined massless fermion perturbed by a field with no static part. 
Although one could eliminate direct corrections to the rest mass by setting S = in 
the hamiltonian, the energies for this case contain terms analogous to the rest energies 
due to a and with a two- fold splitting. The hamitonian is now independent of T3, 
and is invariant under the map 

cr I-)- -cr, T i-)- -T, x H- -x, p t-^ -p. (71) 

The characteristic polynomial of each block Hn is now biquadratic. The four eigenvalues 
are found by choosing one of the alternating signs in the following expression 



E(n) = ±1 



/ b{n) ± a(n) 



(72) 



where 



a{n) = {[l65V(A2a^ + 1) 



+ 



1 + {A + Bya 



{n + 2y- 



8Ba'{A{A + Bya' + A + 2B) {n + 2) (73) 



b{n) 



l + iA' + B'')a' n + 3 1 + iA + Bfa^ 



(74) 



Eigenvectors can be obtained easily. 



6. Entanglement of a confined particle with the external field 

As an example for dynamics in this system we shall now study the evolution of 
entanglement of the Dirac oscillator and the field. We assume our Dirac particle to 
be initially in some eigenstate of the Dirac oscillator, e.g. of positive energy [TT] given 
by 



where 



|x„)=AWH|+) + 4-)|n + l)|-), 



(75) 



n 



n+ {m — 
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e„ = \/2n + l + m? (76) 

The complete system is considered to be initially in a factorized state of the form Dirac 
oscillator x Isospin, i.c-tp = Xn ® X-, X being a Pauli spinor for isospin. We put in 
general |x) = l/v^(cos 6'|+) + sin6'|— )) and examine two different angles 6 = 0,7r/4. 
We shall see that the particular choice of 6 will not be relevant to what we want to 
show. The initial condition can be written as linear combinations of the states (1571) 
with coefficients Ki, i = 1, 4. Such states, in turn, are expressed as superpositions of 
eigenstates of the system \ipD. This results in 



+ i^2E^27'^l^n+2) +^4EM7'^|^n+2), (77) 

where D^^j are coefficients obtained by inverting the expansion 



Applying the spectral decomposition of the propagator to ( 1771) we obtain ip{t). Now 
that we know the explicit form of the time-dependent state, we are in the position to 
compute all quantities related to it. Since we are interested in the way in which the 
bound fermion state entangles with the external degrees of freedom, we take a pure 
state density operator p = \ip{t)){ip{t)\ of the entire system and compute purity P and 
entropy S of the Dirac oscillator subsystem: 



P{t) = TrN,. ((Tr.p(t))') 

S{t) = - TrN,. (Tr,p(t)Log (Tr,p(t))) , 

(79) 

where TrN,(T is the trace with respect to oscillator and *-spin degrees of freedom, while 
Tr,- is the trace with respect to isospin. Results for ( 179|) are depicted in figure ([I]) as 
functions of parameters a, 7, m and time. 

One finds two regimes of maximal entanglement as expected. The first one is 
naturally related to the possibility of tuning the external field to the frequency of the 
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Purity Purity 




Figure 1. Purity and Entropy of the Dirac oscillator state as functions of time t 
(abscissa) and static field 7 (ordinate). Parameters are m = 3.2, a — 1.2. The left 
column corresponds to 9 = tt/A and the right column to 9 — 0. We take n — for 
the Dirac Oscillator state. The structure at 7 ^ m is present in all density plots and 
is persistent at all times. In this region we observe an entropy increase and a purity 
decrease. This is an indication of maximal entanglement between the external field 
and the positive and negative energy states of the confined fermion. The situation 
corresponds to a regime in which the field intensity is comparable to the rest energy: 
particles can be created in the context of quantum field theory. 

oscillator, i.e. a ~ 1. The second regime appears when the strength of the coupling 
is comparable with the rest mass of the Dirac particle, i.e.j ~ m. This shows that 
the structure of the external field is capable of probing the region where particles can 
be created. It is well known that under these circumstances the external field creates 
particles and the quantization of the Dirac field becomes important [21]. At the level of 
classical fields, the results show that this particle creation is accompanied by a surge in 
the entanglement between the Dirac oscillator and the external field. It is not surprising, 
that purity and entropy yield quite similar information. We display both, as entropy 
is more popular but purity can be obtained analytically. We refer to [22], [23] for a 
discussion of the difference of the two quantities. 
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7. Conclusions and outlook 

By considering a bound fermion under the influence of a non-local external fleld, we 
have constructed a generalization of the Dirac oscillator in which interactions represent 
the influence of extended particles in the Yukawa formalism or the non-local gauge field 
theory. The relativistic model has been shown to be solvable in one, two and three 
dimensions by constructing the corresponding invariants in a specific frame of reference. 
By covariance, the solutions in arbitrary inertial frames can be reached by Lorentz 
transformations. We gave spectra and eigenstates explicitly for several cases according 
to different choices of parameters. These include explicit expressions for positive and 
negative energy (particle and antiparticle) states as well as corrections to the rest mass 
of the bound fermion. An interpretation of entanglement of the field with positive and 
negative energy states as a particle creation effect has been possible. 

It is worth to mention that our solvable problem has been formulated in a way 
which allows an immediate identification with quantum optics. The remarkable analogy 
between the Dirac oscillator and the Jaynes-Cummings hamiltonian [17], [18] has been 
pointed out before [24] , [25] , [26] with the aim of producing such a system in a quantum 
optical experiment. Our work may extend the analogy further by identifying our 
hamiltonians with a Jaynes-Cummings model of two atoms. The *— spin and isospin 
of the field correspond to the operators describing two atoms of two levels, while the 
harmonic oscillator operators a, correspond to one mode of the electromagnetic field. 
An experimental realization for this system in the context of quantum optics can be 
proposed as well. In [9], we have given a dynamical study in this context. Specifically, 
we may consider the following setup for the case studied in section 5.1, where the 
parameters are quite general. We require two atoms of different species trapped in an 
electromagnetic cavity. Both atoms must be nearly in resonance with only one mode of 
the eletro magnetic field. Therefore, a sufficiently large spacing between adjacent modes 
is needed. The coupling constant a defined in relation with cu (fixed as unity), can be 
adjusted by placing the atoms in regions with different field intensties. Moreover, a large 
distance between the atoms is needed in order to ignore direct interaction terms. The 
case 5.2 is a particular restriction for which the atoms are of the same kind. Finally, the 
evolution of entanglement studied in section 6 can be realized by preparing the initial 
state in fl77|) . In principle, one could prepare it by trapping one atom in the cavity 
and measuring the energy of the total system (this corresponds to a dressed state of the 
Jaynes-Cummings model with one atom). After this is achieved, the second atom can be 
introduced in the trap. This setup thus emulates our model and makes it experimentally 
accessible. 
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